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- Abstract 

Lev energy proton observations^ ewer the polar cap and in the inter 

planetary and magnetosheath regions are compared during May 26, 1967 * 

The results indicate that at this tine the entry of lov energy solar protons 
into the magnetosphere vas controlled by a diffusion process. It also - 
appears that the protons had access to the entire figure- 8 current pattern 
in the geomagnetic tail and diffused into the two approximately circular 
cylinders containing the field lines frea the northern and southern polar 
caps. A first approximation value of lO^'-lO^cmPsec”* is obtained for 

the' radial diffusion coefficient. D , in the tail. Ko value can he derived 

• . - • w XT 

for the diffusion coefficient along, the field line, D > with the present 

zz 

data. Consequently, in the present case, oily a rough upper l im it of ' 

L s (420-2400)1^ is obtained for the distance to the region where low energy 
protons are diffused into the tail. 



«- 
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We present evidence in this note that diffusion effects can play a 
major role in the entry of lew energy protons into the magnetosphere. 

Tl^se results also indicate that the entire figure-8 pattern of the geo- 
magnetic tail current system f Axford et al . . 1965; Dessler and Juday . 

1965] is rapidly exposed to any available solar proton fluxes. 

The data of interest are simultaneous proton observations in the 
interplanetary medium and over the northern polar cap during the solar 
proton event of May 23, 1967. The interplanetary measurements are from 
the Explorer 3^ satellite and give the intensities of 1-10 Msv protons. 

The polar cap measurements are from satellite 1963 38C and present 
the intensities of 1.2-8. 5 Jfev protons at 1100 km altitude. The local 
time orientation of both satellites is shewn in Figure 1. Data is 
presented only for times when Explorer 3^ was near apogee and statistically 
significant count rates were observed over the polar cap. This period 
is indicated by the solid portion of the Explorer 3^ orbit in Figure 1. 

The time period to be discussed is shown in Figure 2 where both the 
interplanetary and polar cap Froton intensities are plotted. Explorer 34 
data are averaged over 15 minutes. The 13 polar cap passes obtained 
during this period are shown as polar cap averages for invariant latitudes, 
A, greater than J0°. Local times for these passes are also indicated 
and the passes to be discussed in detail are numbered from one to ten. 

For reference a plot of D-- values is included in the figure. 
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In addition we have summarized in Figure 2 the magnetic field con- 
figuration external to the magnetosphere at unis tine. These magnetic 
field data are from Explorer 33 at ~ 0248 hours local tine and at 65 Rg 
and from Explorer 3^ [Fairfield, Behannon, Ness, personal communication; 
Behannon et al . , 1968]. During this time the magnetosphere was imbedded 
in a relatively steady field which changed from strongly southward 
(6 > "45°) to strongly northward (e > 60°) at ~ 1330 hours. 

In spite of an unfortunate gap in the data. Figure 2 indicates that a 
small increase in the intensity of low energy interplanetary protons began 
at ~ 0230 hours on May 2 6. The sequence of five polar cap passes spanning 
this increase (numbers l-y) are shown in Figure 3a. The latitude profiles 
over the polar region are shown along with simultaneous and unnormalized 
interplanetary fluxes. Figure 3b presents similar data for passes six 
through 10. 

Similarities and discrepancies have been observed previously in 
comparing the polar cap and interplanetary time histories of low energy 
proton intensities f Williams and Bostrca . 1967; Krimigis et al . . 

1967]. The absolute unnormalized fluxes shown in Figure 2 point out 
that while general similarities exist between these two regions, significant 
differences in the respective time histories do occur. In particular it is 
seen that, a) the ratio of polar cap to interplanetary fluxes varies by a 
factor of ~ 3 during the period shewn, and b) the significant decrease in 
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the interplanetary fluxes occurring at ~ 1120 hours (Figure 2) is not 
observable over the polar ap except possibly as a slower steady decrease. 
These features do not support the existence at these times of a readily 
accessible polar cap region as in an open magnetospberic configuration. 

They are qualitatively consistent with lew energy proton access to the 
polar cap being controlled by a diffusion process. 

These results are further supported when the individual polar cap 
passes are considered in detail as in Figure 5?- where we show a time sequence 
of five passes (l-5 in Figure 2) obtained during a period in which a slight 
increase occurred in the interplanetary intensities. Here is clearly seen 
a time history of polar cap protons which progresses from pass 1 where the 
intensities tend to increase toward the polar cap, to passes 2, 3> and 4 
which show a profile peaked near the auroral regions and having a minimum 
at the highest latitudes, to pass 5 which shows a distribution much, more 
uniform over all latitudes sampled. Figure 3b continues this polar cap 
time history (passes 6-10) through an additional small interplanetary 
increase occurring at ~ 1520 hours (Figure 2). 

These data are consistent with diffusion into the geomagnetic 
tail being the primary mode of entry of low energy protons into 


the magnetosphere at this time. Such a process has been discussed originally 
by Michel [1965] in considering cosmic ray cutoffs and by Michel and Dessler 


[1965] in considering polar cap absorption inhomogeneities . The general 
features of the polar cap observations can be explained by the magnetotail 


configuration shown in Figure 4a, b. A figure-8 tail current system defines- 
two approximately cylindrical surfaces which when coupled withtthe 
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l rotation of the earth contain the field lines connecting to the lev latitude 

edge of the northern and southern polar caps rDessler and Juday. 1965; 

Axford et al . . 1965]. Diffusion into the tail region should produce polar 
cap profiles which are initially peaked in the auroral regions and which 
gradually fill in the high latitude polar cap region in much the same 
manner as shown in Figure 5. Such profiles are not rare and have 
been reported previously fBostrom et al . , 1967; Zmuda et al . , 1967; Blake 
et al. . 1968]. 

I On the nlghtslde hemisphere the magnetopause field lines for the 

? configuration sketched in Figure 4a connect to the auroral oval. The 

situation is less clear on the dayside hemisphere particularly near the 
noon meridian. During average geomagnetic conditions, field lines from 
the dayside auroral oval do seem to be those passing through the magneto- 
pause neutral region and on back to the tail r Fairfield, 1968]. However, 
it is difficult to move the latitude of the neutral point to the low 
latitudes where aurora and lev energy solar protons are observed during 
magnetically active periods. The exposure of these low latitude 
regions to the tail does seem possible in the nightside hemisphere 
l~ Williams and Hess. 1966; Williams , 1967]. 

A rough estimate of the magnitude of the diffusion coefficient may 
be obtained by considering diffusion into a long cylinder of radius a, 

(Figure 4a), which has an initial internal density of zero ana in which 
diffusion is purely radial. Figure 4c presents the solution to this problem 
by shewing the radial density profile for several different times following 
f the application of an external and steady source of strength n Q f Crank, 1957], 

f ■ 

# 


-5- 


The similarity of the radial profiles in Figure 4c to the polar cap profiles 
of Figure 3 is readily apparent. (Note that r/a = 0 corresponds 
roughly to A 8 */2 81131 r / a * 1 to the nights ide auroral oval. ) The ratio 
of polar cap to auroral oval Intensities will he assumed to yield n/n , 
the ratio of internal to external densities, at r/a = 0, Figure 4c thur. 

p 

yields the appropriate value of Dt/a^ from which D can he obtained. 

Pass 4 in Figure 3 gives a polar cap to auroral region ratio of 
n/n Q «s 0.4. From Figure 4c we obtain a value of Dt/a 2 — 0.15. Due to the 
gap in the Explorer 34 intensities, the appropriate value of t can only 
be estimated from Figure 2 as t as 1.5-5 hours. Using a = 20 Rg gives 
D = 1.4 (10) 15 - 4.0 (10) 15 cm 2 / sec. This range of values represents an 
upper limit to the diffusion coefficient as the polar cap to auroral region 
ratio contains an unknown background intensity from protons present at 
the time of the interplanetary increase. 

Simple diffusion theory gives the following relation .between the diffusion 
coefficient, the mean free path jC, and the particle velocity v. 



If we consider v to be the actual particle velocity, we obtain for a 1 Mev 
proton (v 1.3 (lO)^cm sec -1 ) with the above values of D, 

l - 3(10) 6 - 9(10) 6 cm 

This value is approximately 100 times smaller than the gyroradius, p, of a 

Q 

1 Mev proton in a 20 y field ( p =- 7(10) cm). 
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The above value of t is misleading in the sense that it physically 
applies to radial diffusion (i.e., motion across the field lines) and 
does not take into account the particle’s motion down the field line. 

F -waver, this small value for the mean free path does indicate the 
difficulty of motion across field lines for these low energy protons in 
the tail. This can he seen directly in Vigors 3 where a po].ar 
cap equilibrium distribution is approached only after several hours following 
an interplanetary disturbance. 

An alternative and equivalent description of this diffusion process 
in the tail is obtained by assuming for effective scattering from magnetic 
irregularities* a value of A approximately equal to the gyroradius f Parker , 
19oll. This yields an effective radial velocity, ~ .01 v^ and again 

indicates that these particles travel along the field lines much easier 
than across field lines. 

While these results show the importance of diffusion in the geomagnetic 

tail, a more accurate description requires the use of a tensor rather than 

a scalar diffusion coefficient. However, the principal axes of diffusion 

are defined by the approximate cylindrical geometry of the tail and only 

the components D and D should have to be considered. The present 
rr zz 

15 l6 2 

data indicate that as a first approximation ~ 10 3 - 10 cm /sec 
in the region where diffusion is a major transport mechanism. The relative 
magnitude of D depends on the assumptions employed in estimating pitch 

ZZ 


single scattering effects 



As the above analysis gives no information about D zz , it can give 
no accurate estimate on the "length* of the tail, i.e., the distance to 
the region where particles begin to be diffused into the magnet osphere . 
However, an order of magnitude upper limit may be obtained by assuming 
rectilinear motion of the protons dcwn the "length" of the tail to the 
auroral regions where they first appear. Passes 8 and 10 in Figure 3 
occur after proton intensity increases are observed at Explorer The 

auroral region data of passes 8 and 10 were obtained ~ 10 minutes and 
7 minutes respectively from the attainment of maYirwtm fluxes at 
Explorer 3^» Pass 8 shows a peak intensity in the auroral region and-: thus 
has responded to the interplanetary increase. Pass 10 shows no such 
response. Consequently, assuming rectilinear propagation for 1 Jfev protons, 
we obtain an upper limit for the tail length, L, of 

L s (120 - 2l00)Rg 

Evidence for the existence of a geomagnetic tail has been obtained in 
this range ftfrss et al, 1967 ; Wolfe et al, 1967 ; Marian! and Hess. 1969]* 
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in Figures 2 and. 3 , and reported polar cap electron fluxes rVampola , 

19^91 are also isotropic outside of the loss cone and are indicative of 

diffusive processes. ' •_ 

Delayed access of low energy solar protons to the geomagnetic tail i 

■C"' - r s 

-• ' c : £ 

lias been reported fKane et al, 1968; Montgomery and Singer. 1969]* The ‘ 

lack of a strong energy dependence in the observed delay times [“Montgomery 
and Singer. 1969] o»y he evidence against diffusion effects. However, it 
is not clear that the radial diffusion coefficient should exhibit a 
pronounced energy dependence in the energy range where the proton gyro- 
: radius is much smaller than the radius of the tail (far example, if l 
r ~ constant, D ~ /B ) . 

^ * pT 

Vampola [1969] recently reported observations of energetic 
electrons over the polar regions which . Ja* e characterized by a featureless 
polar cap profile. It is difficult to interpret the electron data, at 
this time in terms of a magnetospberie configuration due to a lack of 
sufficient simultaneous electron data in the interplanetary medium. 

Structural effects might be easily missed if electron diffusion could 
proceed much more rapidly than proton diffusion. Note that the proton 
profiles become' flat and featureless (Passes 5 and 6 , Figure 3) as 

r 

equilibrium is approached. 

The rapid access of low energy protons to the entire figure -8 current i 

£ 

pattern and resulting diffusion into .he cylinders (arrows in Figure lb) can | 

qualitatively account for the app- trance of peak proton intensities in the 
auroral regions and the latitude spread (~ 10°) of this peak. However, it 
should be noted that such a configuration is unable to quantitatively account 
for observed low energy proton cutoffs, particularly on the dayside hemisphere. 

i 
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The ready access of low energy solar protons to 6,6 Rg [Lmaegottl . 
19683 may "be directly related to their appearance in the neutral sheet. 

The diffusion effects discussed by T.angp-nnrf-.fv? [ 1968 ] nay in this case 
be related to diffusion through the "cusp" region and not to radial 
diffusion through the tail. 

In summary we have presented data far a time period in which the 
entry of low energy protons into the geomagnetic tail and their appearance 
over the polar regions was governed by a diffusion process. However, as 
discussed previously rWilliams and Bostrom. 1968], the access of these 
low energy solar protons to the geomagnetic tail and polar cap regions 
depends strongly upon the magnetospheric configuration, which in turn 
■depends upon the external boundary conditions imposed by the interplanetary 
medium. We thus expect to observe strong variations in the manner ir, 
which low energy solar protons enter the magnetosphere . Immediate access, 
delayed access, and even a mixed mode [ “Montgomery and Singer. 19693 are 
possible at different- times. ' 

Of interest is the fact that the present evidence for diffusive 
effects occurred during a period when a strong southward interplanetary 
field was present (Figure 2). It thus seems that a southward interplanetary 
field may not be sufficient cause for immediate access of solar protons 
into the magnetosphere.; 
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Figure Captions 

Explorer 34 and 1963 38 C orbital data for May 25-27, 1967 . 
Solar ecliptic coordinates shown far Explorer 34 orbit. 

Beriod af Interest indicated by heavy line. 

Time history of 1-10 Hev protons observed at high altitudes 
(2 28.2 Eg) and 1.2-8. 5 Mev protons observed at 1100 km 
ever the northern polar cap. The high altitude data are _ 
plotted as 15 -minute averages and the polar cap data are 
averages for invariant latitudes > 7 O 0 . The data shown are 
. absolute and unnormalized fluxes. Orientation of the inter- 
planetary magnetic field with respect to the ecliptic plane 
is indicated. Busfeered polar cap passes are those discussed 
in detail. 

High latitude proton intensity profiles from 1963 3dC for 
(a) passes 1-5 and (b) passes 6-10 from Figure 2. 

Simultaneous Explorer 34 data shown where available. Flux 
values are unnormalized. Explorer 34 altitude > 32 Eg. 

Error bars show statistical uncertainties. 

(a) Sketch of magnet os pher i c configuration; (b) The current 
pattern in the tail region (after Axford et al . , 1965 , and 
Dessler and Juday, 19^5 )• Heavy arrows indicate current 
flow and light arrows indicate proton diffusion; (c) Resulting 
radial diffusion pattern in a long cylinder of radius a, with 
initial internal density of zero, following application of 
an external density, n 0 . The internal density distribution, 
n/n Q , is shown as a function of r/a for various values of 
Dt/a 2 where D = diffusion coefficient and t = time following 
application of the external density, n Q . From (a) it can be 
seen that these radial profiles transform to polar cap 
latitude profiles in the following approximate way: r/a = 0 

corresponds to A = n/2 and r/a = 1 corresponds to the auroral 
region. 
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